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Abstract 

Using examples we test formulae previously conjectured to give the 
fractal information dimension of chaotic repellors and their stable and 
unstable manifolds in "typical" dynamical systems in terms of the Lya- 
punov exponents and the characteristic escape time from the repellor. 
Our main example, a three-dimensional chaotic scattering billiard, yields 
a new structure for its invariant manifolds. This system also provides an 
example of a system which is not typical and illustrates how perturbation 
to the system restores typicality and the applicability of the dimension 
formulae. 

PACS: 05.45. +b 



1 Introduction 

In addition to chaotic attractors, nonattracting chaotic sets (also called chaotic 
saddles or chaotic repellors) are also of great practical importance. In particular, 
such sets arise in the consideration of chaotic scattering, boundaries between 
basins of attraction, and chaotic transients. If a cloud of initial conditions is 
sprinkled in a bounded region including a nonattracting chaotic set, the orbits 
originating at these points eventually leave the vicinity of the set, and there is 
a characteristic escape time, r, such that, at late time, the fraction of the cloud 
still in the region decays exponentially at the rate t _1 . 

The primary focus of this paper will be on studying the fractal dimension 
of nonattracting chaotic sets and their stable and unstable manifolds. Fractal 
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dimension is of basic interest as a means of characterizing the geometric com- 
plexity of chaotic sets. In addition, a knowledge of the fractal dimension can, 
in some situations, provide quantitative information that is of potential prac- 
tical use. For example, in the case of boundaries between different basins, the 
basin boundary is typically the stable manifold of a nonattracting chaotic set, 
and knowledge of the stable manifold's box-counting dimension (also called its 
capacity) quantifies the degree to which uncertainties in initial conditions result 
in errors in predicting the type of long-term motion that results (e.g., which 
attractor is approached; see Sec. 5 and Ref. Q). Our focus is on obtaining the 
information dimension of a suitable "natural measure" fi lying on the chaotic 
set (see Sec. 2 for definitions of the natural measures for nonattracting chaotic 
sets and their stable and unstable manifolds). The information dimension is a 
member of a one parameter (q) class of dimension definitions given by 

Dq = ]im J- l ^d, (i) 

q e^ol-q ln(l/e) V ; 

where q is a real index, e is the grid spacing for a (i-dimensional rectangular 
grid dividing the d-dimensional state space of the system, and /ii is the natural 
measure of the grid cube. The box-counting dimension is given by ([!]) with 
q = 0, and the information dimension is given by taking the limit (7 > 1 in (|lj) , 

Dl = 1^0 hjlje) ' /(e) = ^ Mi ln(1/Vi) ' (2) 

In general, the information dimension is a lower bound on the box-counting di- 
mension, Di < Dq. In practice, in cases where D\ and Dq have been determined 
for chaotic sets, it is often found that their values are very close. 

We are specifically concerned with investigating formulae conjectured in Ref. 
D that give the information dimensions for the nonattracting chaotic set and its 
stable and unstable manifolds in terms of Lyapunov exponents and the decay 
time, t. These formulae generalize previous results for nonattracting chaotic 
sets of two-dimensional maps with one positive and one negative Lyapunov 
exponent H || , and for Hamiltonian systems of arbitrary dimensionality (s| . In 
turn, these past results for nonattracting chaotic sets were motivated by the 
Kaplan- Yorke conjecture which gives the information dimension of a chaotic 
attractor in terms of its Lyapunov exponents |J. A rigorous result for the 
information dimension of an ergodic invariant chaotic set of a two-dimensional 
diffeomorphism has been given by L.-S. Young Q, and this result supports the 
Kaplan- Yorke conjecture for attractors and the two-dimensional map results of 
Refs. H and Q for nonattracting chaotic sets. Furthermore, Ledrappier j|] has 
proven that the conjectured Kaplan- Yorke dimension formula is an upper bound 
on the information dimension of chaotic attractors. Soluble examples for typical 
systems and numerical experiments tend to support the Kaplan- Yorke result. 
On the other hand, some specific examples violating the Kaplan- Yorke formula 
also exist. Thus, in the Kaplan- Yorke conjecture they claim that their formula 
for the dimension applies for "typical" systems. That is, if a specific system 
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violates the formula then by making a "typical" arbitrarily small change to the 
system we create a new system for which the Kaplan- Yorke formula holds (this 
is discussed in Refs. [|| and p0|). Specific systems violating the conjecture are 
called atypical, and, by Ledrappier's result, have dimensions below the Kaplan- 
Yorke prediction. The conjectured results of Ref. for nonattracting chaotic 
sets are also of this character in that they are also claimed to hold only for 
"typical" systems. Thus, the issue of what constitutes a typical system (or a 
typical perturbation of an atypical system) is central. At present there is no 
rigorous formulation of "typical" for this purpose. For this reason it is important 
to address this question through examples. One of the purposes of this paper is 
to do just that. More generally, although many examples, both analytical and 
numerical, exist which support the dimension formula for chaotic attractors, for 
the case of nonattracting chaotic sets there are virtually no supporting examples 
(except for the case of two-dimensional maps with one positive and one negative 
Lyapunov exponent ||[ |J). This paper provides such examples. 

In Section 2 we review the results of Ref. [0. We define the natural tran- 
sient measure on nonattracting chaotic sets and on their stable and unstable 
manifolds, and give the conjectured information dimension formulae for these 
measures. 

In Section 3 wc introduce a noninvertible two-dimensional map with two 
expanding (positive) Lyapunov exponents. This map is simple enough that it 
can be fully analyzed. This analysis is especially useful in elucidating certain 
aspects of the general problem, particularly in relation to the natural transient 
measure and the role of fluctuation in the finite time Lyapunov exponents. 

In Section 4 we introduce our main example, a hard walled billiard chaotic 
scattering system in three spatial dimensions. For a particular symmetric con- 
figuration of the scatterer, we find (Sees. 5 and 6) that the situation is atyp- 
ical (i.e., the conjectured formula does not apply), but that deviation of the 
configuration from the symmetric case immediately restores the validity of the 
conjectured dimension formula. 

Section 5 presents numerical results for the dimension in the atypical (sym- 
metric) and typical (asymmetric) cases obtaining good agreement with theoret- 
ical results in both cases. 

Section 6 gives a theory for the structure of the nonattracting chaotic in- 
variant set present in examples modeling our billiard scatterer. It is shown that 
the stable manifold of the invariant set is a continuous, nowhcrc-diffcrcntiable 
surface in both the typical and the atypical cases. Furthermore, we are able to 
obtain an explicit formula for the dimension in the atypical case (this formula 
is compared to numerical computations in Sec. 5). 

In Sec. 7.1 we present a derivation of the dimension formulae of Sec. 2. 
Although these formulae have been previously derived in Ref. Q , the previous 
derivation assumed a specific type of mapping. The situation envisioned in Ref. 
[^| does not apply to our billiard example. Thus, we are motivated to provide 
more general arguments for the applicability of the formulae in Sec. 2. In 
Sec. 7.2 we derive a formula for the dimension of the stable manifold of the 
chaotic repcllor found in the typical (asymmetric) case. This derivation shows 



3 



explicitly where the derivation of the typical formula fails for this nontrival 
atypical system. 

In Sec. 8 we present conclusions and further discussion. In particular, we 
point out that the continuous, nowhere-differentiable surface forming the stable 
manifold of our billiard example may be amenable to experimental observation 
in a situation where the billiard walls are mirrors and the orbits are light rays. 
In general, we comment that we believe that such optical billiard experiments 
provide a particularly convenient arena for the experimental investigation of 
different types of fractal basin boundaries. This has so far proven to be very 
difficult in other experimental settings. 



2 Dimension Formulae 

A chaotic saddle, A, is a nonattracting, ergodic, invariant set. By invariant we 
mean that all forward and reverse time evolutions of points in A are also in A. 
The stable manifold of A is the set of all initial conditions which converge to 
A upon forward time evolution. The unstable manifold of A is the set of all 
initial conditions which converge to A upon reverse time evolution. We say A is 
nonattracting if it does not completely contain its unstable manifold. In such a 
case there are points not in A that converge to it on backwards iteration. 

To define the characteristic escape time, r, first define a bounded region, 
R, which contains A and no other chaotic saddle. Uniformly sprinkle a large 
number, N(0), of initial conditions in R. (In this section we take the dynamical 
system to be a discrete time system, i.e., a map.) Iterate the sprinkled initial 
conditions forward n» 1 times and discard all orbits which are no longer in R. 
Denote the remaining number of orbits N(n). We define r as 



N(0) 



(3) 



or, more formally, r = lim^oo limjv(o)^oo hi[N(0)/N(n)]/n. The Lyapunov ex- 
ponents are defined with respect to the natural transient measure of the chaotic 
saddle JllJ . This measure is defined on an open set C C R as 

KC) = lim hm N{ [ n > n : C \ (4) 

V ' n^oojv(0)^oo N(n) 

where < £ < 1, and N(m,n,C) is the number of sprinkled orbits still in R 
at time n that are also in C at the earlier time m < n. The above definition 
of /x(C) is presumed to be independent of the choice of £ as long as < £ < 1 
(e.g., £ = 1/2 will do). 

We take the system to be M-dimensional with U positive and S negative 
Lyapunov exponents measured with respect to /x (where U + S = M) which we 
label according to the convention, 

h% > h^ J _ 1 > ■ ■ ■ > h+ > > -hi >■■■> -hg_ x > -hg. 
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Following B] we define a forward entropy, 

U 



We now define a natural transient measure /is on the stable manifold and a 
natural transient measure fj,u on the unstable manifold. Using the notation of 
Eq- (@), 

, S (C) = lim hm (5) 

V ' "^coAT(o)^oo N(n) y ' 

fiu(C) = lim lim N{l ?;J l \ C) ■ (6) 

^ V ^ n->oo JV(0)->oo iV(n) 

Thus, considering the N(n) orbits that remain in i? up to time n, the fraction 
of those orbits that initially started in C gives fis(C), and the fraction that end 
up in C at the final time n gives fj,u(C). We use the measure (^), (|^), (|J) to 
define the information dimensions of the invariant set, the stable manifold, and 
the unstable manifold, respectively. 

According to Ref . M , the dimension of the unstable manifold is then 

D u = U + I + H -« l \ + - + h '\ (7) 
h i+i 

where I is defined by 

H h hj + hj +1 >H>h^-\ +hj. 

The dimension of the stable manifold is 

IW + J+ '-Wf"^*) , (8) 

where J is defined by 

/i+ + • • • + h+ + h+ +1 > H > h+ + ■ ■ ■ + h+. 

Considering the chaotic saddle to be the (generic) intersection of its stable and 
unstable manifolds, the generic intersection formula gives the dimension of the 
saddle, 

D A = D v + D s - M. (9) 

It is of interest to discuss some special cases of Eqs. (|7|)-@. In the case of a 
chaotic attractor, the invariant set is the attractor itself, the stable manifold is 
the basin of attraction, and we identify the unstable manifold with the attractor. 
Thus Ds = M and D\ = D\j. Since points near the attractor never leave, we 
have t — co. Equation (0) then yields the Kaplan- Yorke formula ||, 

(ht H h hi T ) - (hT + . . . + h7) , , 

D A = U + I+^ u> _ \ i 1_±^ ( 1Q ) 

h i+i 
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where I is the largest integer for which (hf + . . . + h^) — (hj + . . . + hj) is 
positive. 

In the case of a two-dimensional map with one positive Lyapunov exponent 
ht and one negative Lyapunov exponent h± with the exponents satisfying ht — 
K - 1/r < 0, Eqs. (0) and (g) give the result of Ref. § and gj, 



D rj = 1 



Do = 1 



hf - 1/r 



h i 



h+ - 1/r 



Another case is that of a nonattracting chaotic invariant set of a one-dimensional 
map. In this case 5 = and {7 = 1. The unstable manifold of the invariant 
set has dimension one, Dy — 1. Recalling the definition of the stable manifold 
as the set of points that approach the invariant set as time increases, we can 
identify the stable manifold with the invariant set itself. This is because points 
in the neighborhood of the invariant set are repelled by it unless they lie precisely 
on the invariant set. Thus, D$ = -Da, and from (||) and (^|) we have D$ = D\ = 
H/ht, where H = ht - 1/r. 

Still another simple situation is the case of a two-dimensional map with two 
positive Lyapunov exponents. This case is particularly interesting because we 
will be able to use it (Sec. 3) to gain understanding of the nature of the natural 
measure whose dimension we are calculating. In this case U — 2 and S = 0. 
Thus Djj — 2 and Dg = D^. There are two cases [corresponding to J = and 
J = 1 in Eq. (0)]. For K^t < 1, we have that D$ = -Da is between zero and 
one, 

Ds = D A = l + ^--^. (11) 
For h^r > 1, we have that D$ = -Da is between one and two, 

D s = D A = 2--±-. (12) 

hjT 

In the next section we will be concerned with testing and illustrating Eqs. 
and (O) by use of a simple model. 



3 Illustrative Expanding Two- Dimensional Map 
Model 

We consider the following example, 

x n +i = 2x„ modulo 1, (13) 
Vn+l = K x n)Vn + tt - sin(27T2; n ), (14) 



G 



where X(x) > 1, and the map is defined on the cylinder — oo < y < +00, 
1 > x > 0, with x regarded as angle-like. We take A (a;) to be the piecewise 
constant function, 

A(x) = { ^ °<* <1/2 : (15) 
w 1 A 2 1/2 < x < 1, K J 

and, without loss of generality, we assume Ai < A 2 . [Later, in Sec. 3.5, we will 
consider the problem with a general function X(x) and a general chaotic map 
x n +i = M(x n ) replacing (|l3|), but for now we focus on (|l3|)-(p"5|).] 

For this map, almost every initial condition generates an orbit that either 
tends toward y = +00 or toward y = — 00. Figure |l| shows the regions where 
initial conditions generate these two outcomes, with the black (white) region 
corresponding to orbits that tend toward y — —00 (y = +00). Initial conditions 
on the border of these two regions stay on the border forever. Thus, the border 
is an invariant set. It is also ergodic by virtue of the ergodicity of the map 



x n+ i = 2x n modl. We wish to apply Eqs. (11) and ( |12[ ) to this invariant set 
and its natural measure. 

The Jacobian matrix for our model is 



J(x) 



2 

77COs27ra; A(x) 



Thus, for an ergodic invariant measure of the map, the two Lyapunov exponents 
are 

h a =plnAi + (1 -p)lnA 2 (16) 

and 

h b = In 2, 

where p is the measure of the region x < 1/2. To find h a we thus need to know 
the measure of the invariant set. The measure we are concerned with is the 
natural transient measure introduced in Sec. 2. 



3.1 The Decay Time and the Natural Measure 

Consider a vertical line segment of length Iq whose x coordinate is xo and 
whose center is at y = yo- After one iterate of the map (p^|)-jl5|), this line 
segment will have length l\ = A(xo)4> and be located at x = X\ with its center 
at y = j/i, where (x\,y\) are the iterates of (£0,2/0) using the map (|l3|)-(|l5|). 
Thus we see that vertical line segments are expanded by the multiplicative 
factor \(x) > Ai > 1. Now consider the strip, —K < y < K, and sprinkle 
many initial conditions uniformly in this region with density p . A vertical line 
segment, x — xq, —K < y < K , iterates to x = X\ and with its center at 
2/i = {f]/2n) sin27ra:o. We choose K > (ry/27r)(Ai — l)" 1 so that the iterated line 
segment spans the strip —K<y<K. After one iterate, the density will still be 
uniform in the strip: The region x < 1/2 [x > 1/2), —K < y < K 7 is expanded 
uniformly vertically by Ai (A2) and horizontally by 2. Thus, after one iterate, 
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the new density in the strip is p\ = [(X 1 1 + A 2 1 )/2]po, and, after n iterates, we 
have 

p„ = [(Ar 1 +A2 1 )/2] n po- 

Hence the exponential decay time for the number of orbits remaining in the 
strip is 

1 



In 









[i 


VAi " 





(17) 



(n) 

To find the natural stable manifold transient measure of any x- interval s m — 
[m/2™, (m+l)/2 n ], where m = 0, 1, . . . , 2™ — 1, we ask what fraction of the orbits 
that were originally sprinkled in the strip and are still in the strip at time n 
started in this interval. Let Sm experience n\(m) vertical stretches by Ai and 
ri2(m) = n — ni(m) vertical stretches by A2. Then the initial subregion of 
the st' still in the strip after n iterates has vertical height KX[ nit " m ^ X^ 112 ^ ■ 
Hence the natural measure of st' is 

— n \ — ni(m) ■, -n 2 (m) ^n 2 {m) , ni(m) 

U (M)) - 1 A i 21 - 21 21 Ma) 

M m ] ~ [K A r x +A2 1 )]™ ^ (Ai + a 2 )« • 

(Note that this is consistent with fi([Q, 1]) = XH=oM s ™') = Thus the 
measures of the intervals [0, 1/2] and [1/2, 1] are 

and 

1- p = a*(*i ) = i—^r- 

M + A2 

It is important to note that our natural transient measures p and (1 — p) for 
the two-dimensional map are different from the natural measures of the same 
x-intervals for the onc-dimcnsional map, x n +i = 2x n mod 1, alone. In that case, 
with probability one, a random choice of x$ produces an orbit which spends half 
its time in [0, 1/2] and half its time in [1/2, 1], so that in this case the natural 
measures of these regions are p = (1 — p) = 1/2. The addition of the y-dynamics 
changes the natural measure of x-intervals. 

From (|l6|) we obtain 

A - ■ In Ai + Al InAa. (19) 



Ai + A2 Ai + A2 

For a general function f(Z) with cPf/dZ 2 < 0, averaging over different values 
of Z gives the well-known inequality (f(Z)) < f((Z)) where ((• • •)) denotes the 
average of the quantity (• • •). Using f(Z) —hiZ with Z = Ai with probability 
p = A2/ (Ai + A2) and Z — A2 with probability (1 — p), this inequality and Eqs. 
(0) and |H) yield the result that 

K<-. (20) 
r 
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In fact, we will see in Sec. 3.5 that (|2(j) remains true for any choice of the func- 
tion X(x) > 1 in (|l4|). In ( p(i| ) the equality sign applies if the vertical stretch- 
ing is uniform (Ai = A2) but, for any nonuniformity in the vertical stretching 
(Ai 7^ A2), h a is strictly less than 1/t. 



3.2 Application of the Dimension Formulae 

Let A2 = rAi, r > 1, and imagine that we fix r and vary Ai. Applying Eqs. ([11 
and ( |l2| ) to our example we obtain three cases, 

(a) h b > 1/t > h a (Ai small), 

(b) 1/t > h b > h a (Xi moderate), and 

(c) 1/r > h a > h b (Ax large). 



Corresponding to these three cases ( pd[ ) and (12) yield the following values for 
Da, the dimension of the invariant set, 

lnfl + r" 1 ) -InAi , x 

D a = l + — j-^ -,forAi<A„, (21) 



m Ai + (1 + r) 1 Inr 



ln(l + r- r ) + (l + rj-Mn 



r 



D c = ^ ^ '- ,fijrA 6 <Ai ) (23) 

where In A a = ln(l + r _1 ) and In Af, = In 2 — (1 + r)^ 1 lnr. The solid line in Fig. 
H shows a plot of Da versus In Ai for r — 3. Note that for large Ai, Da — D c is 
independent of Ai. 

It is also instructive to consider the case of uniform stretching (r = 1) for 
which Ai = A2. In that case, h a = 1/t, and there is a rigorous known result for 
the dimension 0. For Ai = A 2 , Eqs. ©-(H) yield 

DA -f 2 -T^ fOT 1<A X <2, 

y 1 for Ax > 2. [M) 

(For r — * 1 region (b), where Da = D b , shrinks to zero width in Ai.) For 
r = 1 the natural transient measure is uniform; from ( |l8|) we have fi(sm ) — 
2~ n independent of the interval (i.e, independent of m). In this case there 
is no difference between the capacity dimension of the invariant set and the 
information dimension of its measure. Equation (^4|) agrees with the rigorous 
known result, thus lending support to the original conjecture. 



3.3 Numerical Tests 

The formulae (pl|)-(p3"|) were verified by numerical measurements of the infor- 
mation dimension, D\, of A at various values of Ai with r = X2/X1 fixed at 
r = 3 (Fig. |^). Shown for comparison is the box-counting dimension, Dq. The 
values of the box-counting dimension are numerically indistinguishable from the 
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values of the information dimension when Dq, D\ > 1, or Ai < X a , the region 
corresponding to formula (f2l|). For Ai > A a , A is a smooth curve and so has a 
box-counting dimension of Dq — 1. No points for D\ are shown near Ai = A&. 
It can be argued (Appendix A) that numerical convergence is too slow here to 
yield accurate measurements of the dimension. 

To numerically determine the information dimension of A (the data shown 
as open squared in Fig. ^) , we place a square x — y grid with a spacing e between 
grid points over a region containing A. Using the method dscribed in the next 
paragraph we compute the natural measure in each grid box and repeat for 
various e. The information dimension is then given by 

A = lim 



♦o ln(l/e)' 

where /(e) = J^iLi l n (V Mi) is a sum over the N(e) grid boxes which intersect 
A and fii is the natural measure in the ith box. The slope of a plot of 1(e) versus 
In e gives D\ . The box-counting dimension (the data shown as black dots in Fig. 
||) is given by 

-Do = hm 

e^o lnl/e 

and calculated in an analogous way. 

To determine which boxes intersect A and what measure is contained in each 
of them we take advantage of the fact that A is a function 111]] [see Appendix 



A, Eq. (A2)]. That is, for each value of x there is only one corresponding value 
of y in A, which we denote y = y\(x). We divide the interval < x < 1 into 2 n 
intervals of width S = 2~ n . We wish to approximate y\(xo) for in the center 



of the a;- interval. To do this we iterate xq forward using (13) m times until the 
condition 

6 -XT l XT > 1 (25) 

is first met, where mi (777-2) is the number of times the orbit lands in0<a;<l/2 
(1/2 < x < 1), and mi + m 2 = m (we will see the reason for this condition 
below). All of the values, Xi, of the iterates are saved. Starting now from x m 
and taking y m = we iterate backward m times. For 77 small enough, A is 
contained in — 1 < y < 1, so the point (x m , y m — 0) is within a distance 1 in the 
y-direction of A. The m reverse iterations shrink the segment y m < y < yA(x m ) 
by a factor A™ 11 A™ 2 so that \yo — y\(xo)\ < S/2 by condition (p5|). Thus, we 
have found a point yo that approximates y\(xo) to within S. Since we will 
be using e boxes with e ^ 5, we may regard yo as being essentially equal to 
Va(%o)- The measure in the 6 width interval containing a;o is found by iterating 

Xq forward 77 times and using equation (|1§|), /i = fX[+%fj" > wnere n i ("2) is the 
number of times the orbit lands in < x < 1/2 (1/2 < x < 1), and n\ + 112 = n. 
We associate this measure with the point (2:0,7/0) (with yo found by the above 
procedure) . 
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Note that for fractal y\(x) the y interval occupied by the curve y = y\(x) 
in an x interval of width S <C 1 is of order S 00-1 which is large compared to 
5. We now cover the region with new grids having successively larger spacing, 
£, = 2 l 5 = 2 l ~™, and calculate I(si) and N(si) based on the data taken from 
the first 5-grid. For i large enough, such that the y extent of the curve y\(x) 
in a typical S width interval is less than e,; (i.e., ej > 5 2 ~ D ° or i > (Do — l)n) 
we observe linear scaling of log/(ej) and logN(ei) with loge,, and we use the 
slope of such plots to determine D\ and Dq. The dimensions D\ and Dq are 
then determined as described above. 



3.4 Atypical Case 

The conjecture of Ref. is that the dimension formulae of Sec. 2 apply for 
"typical" systems. To see the need for this restriction consider Eqs. (|l4|) for 
the case where 77 = 0. It is easily shown that the dimension formulae can be 
violated in this case. The claim, however, is that rj = is special, or "atypical" , 
in that, as soon as we give 77 any nonzero value, the validity of the dimension 
formulae is restored. In this connection it is important to note that as long as 
r\ 7^ 0, the dimension of the invariant set is independent of the value of 77. This 
follows since if 77 7^ we can always rescale the value of 77 to one by the change 
of variables y = yjr\. To see the violation of the dimension formulae for 77 = 0, 
we note that in this case, by virtue of jl4|), the line y = is invariant. Thus 
the measure is distributed on a one-dimensional subspace, the ic-axis. (This is 
very different from the picture in Fig. |l|, where the invariant set, the boundary 
between black and white, appears to be fractal.) Using the definition of the 
information dimension and dividing the x-axis into intervals of width 2 _n , the 
information dimension of the natural measure is 

Da _ lim ss^)Mwa. (26) 

«^oo ln(2") v ; 



The quantity whose limit is taken in (|26|) is in fact independent of n. Thus, 
taking n = lwe obtain for Da the result that, for rj = 0, 

D A = D c , 

for all Ai and A2 > 1, where D c is given by (^3|). Thus, for h a < hb, Da is 
greater when 77 7^ than when 7/ = 0, and, thus, the conjectured stable manifold 
dimension formula of Sec. 2 is violated. For h a > hb, Da is the same in both 
cases. 

In Sec. 4 we consider chaotic scattering in a three-dimensional billiard ex- 
ample for which the character of the atypical case is more interesting than in 
the above example. In particular, we find for our billiard that in both the typ- 
ical and the atypical cases the stable manifold can have noninteger capacity 
dimension. 
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3.5 General Considerations 



The previous considerations readily generalize to the case of an arbitrary smooth 
function X(x) > 1 and a general chaotic map, x n+ \ — M(x n ), which replaces 
(3.1). Consider the finite time vertical Lyapunov exponent, 

1 n 

h(x,n) = -J2 lnA(M m - 1 (x)) 

m— 1 

computed for the initial condition x. Choosing x randomly with uniform prob- 
ability distribution in the relevant basin for chaotic motion [e.g., x in [0, 1] for 
Eq. (|l3|)], h(x, n) can be regarded as a random variable. Let P(h, n) denote its 
probability distribution function. For large n, we invoke large deviation theory 
to write P(h, n) as |ll| 

\nP(h,n) = -nG{h) + o{n), 

or, more informally, 

P(h,n) ~ e - nGW , (27) 

where the specific form of G(h) depends on M(x) and the specific X(x), and 
G(h) is convex, d 2 G(h)/dh 2 > 0. For the normalization, J P(h,n)dh = 1, to 
hold for n — > oo, we have that 

min G(h) = 0; 

h 

see Fig. |^, where h denotes the value of h for which the above minimum is 
attained. As n — > oo we see that P approaches a delta function, 5{h — h). Thus, 
h is the usual infinite time Lyapunov exponent for almost all initial conditions 
with respect to Lebesgue measure in < x < 1. 

As described above P(h, n) is the probability distribution of h{x 1 n) for x 
chosen randomly with respect to a uniform distribution in [0, 1]. We now ask 
what the probability distribution of h(x, n) is for x chosen randomly with respect 
to the natural transient measure for our expanding map, x n +\ = M(x n ) and 
(pT|). To answer this question we proceed as before and consider an initial 
vertical line segment | gy | < K starting at x [with K > (?7/27r)(A m i n — 1) _1 , 
A m in = muxr X(x) > 1]. After n iterations, this line segment lengthens by the 
factor exp[nh(x, n)J . Thus, the fraction of the line still remaining in the strip 
\y\ < K is exp[— nh{x, n)]. Hence, the fraction of points sprinkled uniformly in 
the strip that still remains after n iterates is 

e -n/T _ J e -nG(h)-nh dhj ( 2 g) 

and the probability distribution of finite time vertical Lyapunov exponents for 
x chosen randomly with respect to the natural transient measure is 

g — nG(h) — nh 

fe -nc ih) - nhdh - (29) 
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Evaluating @ for large n we have J e^^^d/i - e -n[G(/ l ,)+/ l .] ) where 
min[G(/i) + /i] = G(h*) + h* and /i* is the solution of dG(h*)/dh* = — 1. Thus . 



1/r = G(/i*) + /i*. 



(30) 



See the construction in Fig. [| in which the dotted line of slope — 1 is tangent 
to the graph of G(h) at the point h — /i*. The infinite time vertical Lyapunov 
exponent for the transient natural measure is 



h a 



hP(h,n)dh. 



(31) 



Using (HI) and again letting n be large ([II]) yields h a — h*. Referring to Fig. || 
we see that 

h a < 1/t. 

That is, Eq. (£0) is valid for general M(x) and X(x) and not just for M(x) and 
A(a;) given by (|13|) and (15). 



4 A Three-Dimensional Billiard Chaotic Scat- 
terer 

We consider a three degree-of-freedom billiard. The billiard (Fig. ||) is formed 
by a hard ellipsoid of revolution, placed in a hard, infinitely long tube with cross- 
section as shown in Fig. ||(b). The center of the ellipsoid is placed at the center 
of the tube. We consider two cases: (a) the major axis of the ellipsoid coincides 
with the z-axis [see Fig 0(a)], (b) the major axis of the ellipsoid lies in the y-z 
plane and makes an angle £ with the z-axis. The ratio of the minor radius of 
the ellipsoid (rii ) to the width of a side of the tube [dashed line in Fig. |i](b)] 
is 1/4. This leaves the major radius, rj_, and the tilt angle, £, as parameters. 
A point particle injected into the system experiences specular reflection from 
the ellipsoid and the walls (i.e., the angle of reflection is equal to the angle of 
incidence, where both are taken with respect to the normal to the surface off of 
which the particle bounces). When the orbit has passed the top (bottom) of the 
ellipsoid, with positive (negative) z-velocity, we say that it has exited upward 
(downward). We fix the conserved energy so that \v\ = 1. 



4.1 Pictures of the Stable Manifold 

We claim that the case of the untiltcd ellipsoid [£ = and case (a) above] is 
atypical in the sense that it violates the formulae of Sec. 2. However, as soon 
as £ 7^ 0, we claim that the formulae of Sec. 2 apply. We begin by discussing 
the atypical (untilted) case. 

By the symmetry of the geometry of the billiard with the ellipsoid axis 
along z [Fig. f|(a)], the chaotic saddle, A, of this system is the collection of 
initial conditions satisfying z = v z =0. Started with these initial conditions, an 



13 



orbit will have z = v z = for all forward and reverse time. The surface normals 
of the ellipsoid and walls at z — lie in the z = plane, and thus the particle 
cannot acquire a non-zero v z . The z = slice through the three-dimensional 
billiard is a two-dimensional billiard with concave walls [Fig. [|(b)]. It is known 
[pi that a typical orbit in this billiard will fill the phase space ergodically. Near 
z = 0, we can picture a typical point on the stable manifold (denoted SM) 
as having, for example, z slightly less than zero and v z slightly greater than 
zero. The particle will hit the ellipsoid below its equator and, thus, v z will be 
decreased with each bounce, yet remain positive. With successive bounces, the 
orbit on SM slowly approaches z = v z = 0, the chaotic saddle. 

To visualize SM, we note that it forms the boundary between initial con- 
ditions which escape upward and those which escape downward. We say that 
points which, when iterated, eventually escape upward (downward) are in the 
basin of upward (downward) escape. Points which are on the boundary be- 
tween the two basins never escape at all, i.e. they are in SM. We initiate a 
(two-dimensional) grid of orbits (500x500) on the plane,— 3 < x < 3, y = 5.1, 
—2.5 < z < 0, v x = 0, v z = .1, and v y is given by the condition \v\ — 1. We 
iterate each of these initial conditions forward until it escapes. Then we plot a 
white (black) point for each orbit in the upward (downward) basin. The result is 
shown in Fig. ||(a) ■ The boundary between the white and black regions is then 
the intersection of SM lying in the phase space of the five-dimensional billiard 
(x, y, z, v Xl v y , v z constrained by |w| = 1) with the specified two-dimensional x, z- 
plane. SM appears to take the form of a nowhere-differentiable curve. This is 
true in various 2D slices, none of which are chosen specially, which suggests that 
SM has this form in a typical slice. A similar procedure can be followed for the 
case of the tilted ellipsoid and the resulting picture is shown in Fig. ||(b) for 
the case of a tilt angle of £ = 27r/100 radians. 

4.2 Lyapunov Exponents, Decay Times, and Approximate 
Formulae for the Stable Manifold 

Again, we begin with the untiltcd case. To construct a map from this system 
we record the cylindrical coordinates (z, 4>) and their corresponding z and <j) 
velocity components, which we denote (v, lo), each time the particle hits the 
ellipsoid. The coordinate r is constrained, for a given z, by the shape of the 
ellipsoid surface, and v r is given by the energy conservation condition \v\ = 1. 
The four components (z, v, <f>, uj) n give the state of the system at discrete time 

z 

n, where n labels the number of bounces from the ellipsoid. Let z = 



and <j) = 
notation, 



For future reference, we express the map using the following 




M z (z n ,4> n ), 

M^Zn, (j> n ). 



(32) 
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In what follows, when we refer to an orbit, saddle, invariant set, stable manifold, 
etc., we are referring to these quantities for the discrete time map (rather than 
the original continuous time system). 

In the case of the untilted ellipsoid, linearizing about an orbit on A, (i.e., 
z n = 0, 4> n )i we obtain, for the evolution of differential orbit perturbations 5z 
and 6(f), 

Sz n+1 = DM z (0,(j) n )Sz n , 

where $ n +i = M^(0, $ n ) is the map for the two-dimensional billiard of Fig. |i](b), 
DM z (0, cj)) is the tangent map for differential orbit perturbations in z evaluated 
at z = 0, and DM^lO, </>) is the tangent map for differential perturbations lying 
in A. Let ±/i z and ±h<p denote the Lyapunov exponents with respect to the 
natural transient measure (Sec. 2) for perturbations in z and in 0, respectively 
(these exponents occur in positive-negative pairs due to the Hamiltonian nature 
of the problem) . 

In principle, one could numerically evaluate h z and by sprinkling a large 
number, N, of initial conditions in the vicinity of A, iterating n 3> 1 times, 
evaluating h z and over those orbits still near A, and averaging h z and h$ 
over those orbits. We could also find r by this procedure; it is the exponential 
rate of decay of the orbits from the vicinity of A. For cases where the escape 
time r is not long, this procedure, however, becomes problematic. For finite N 
the number of retained orbits can be small or zero if n is too large. Thus, we 
adopt an alternate procedure which we found to be less numerically demanding. 

In particular we make user of the ideas presented in Sec. 3.5. We define the 
uniform measure as the measure generated by uniformly sprinkling many initial 
conditions in A (the hyperplane z = v z = 0). An average over these orbits of 
the tangent space stretching exponents would yield uniform measure Lyapunov 
exponents. We denote the distribution of finite-time Lyapunov exponents with 
respect to this measure by P(h^,, h z ,n) ~ e -nG(h$,h*) j as m (g^)] ) where the 
tilde indicates finite time exponents for initial conditions distributed according 
to the uniform measure. 

To compute the decay time and the Lyapunov exponents with respect to 
the natural transient measure we note that orbits near a point <f> in A iterate 
away from A as exp[n/i z (0, n)]. Thus, the fraction of a large number of initial 
conditions sprinkled near A which remain near A after n iterates is 

J P(J H ,h z ,n)e- nh *dh z , 

and htf, (the infinite time Lypunov exponent with respect to the natural transient 
measure) is 

Jh 4> PCh^,h z ,n)e' n '^dh 
hi, = hm = — = = . 

n^co J P(h 4n h z ,n)e~ nh -dh 
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This expression can be approximated numerically by choosing N initial condi- 
tions, (pi, uniformly in A and calculating 



nh. 



We calculate the finite-time Lyapunov exponents h^((j)i,n) and h z (cj)i,n) using 

the QR decomposition method . Since we chose the <j>i uniformly in A these 

exponents are distributed according to P(n, h^, h z ). N is taken to be large 

enough that there are at least 100 terms contributing to 90% of each sum. The 

range of n is from 10 to about 40. We find that ( nh^ \ versus n is well-fitted 

. n 

by a straight line and we take its slope as our estimate of h^. 

To find r, first note that we numerically find that does not vary much 
from h$ (.91 < h^/h^ < 1) as we change the system parameter (the height 
of the ellipsoid), where the overbar denotes an infinite-time Lypunov exponent 
with respect to the uniform measure on A. Thus, we make the approximation 
P(h ( f,,h Zl n) w P(h z ,n) 7 G(h^,,h z ) w G(h z ) and plot G(h z ) versus h z . The 
value of 1/r is given by Eq. (|30|), and h z is given by dG(n,h z )/dh z = — 1. 
(Note that h z < 1/t < h z .) A third order polynomial is fit to the data for 
G(h z ) and used to find r and h z . This calculation is performed at a value of n 
which allows a significant number of points in G(h z ) versus h z to be collected 
near in the range h z < h z < h z . 

For the case of the tilted ellipsoid, we will consider a very small tilt angle, 
£ = 27r/100. With this small tilt angle, the Lyapunov exponents and decay 
times for the tilted and the untilted cases are approximately the same. Thus, 
for the tilted case, we will use the same Lyapunov exponent values, ±h z and 
±/i0 and decay time r, that we numerically calculated for the untilted case. The 
dimension formula of Sec. 2 for SM is 

Ds = 4- {h^Y 1 for/i0T>l, , . 

D s <i farV<l. 

(As for cases (b) and (c) of Sec. 3.2, we could write explicit expressions for Ds 
for hff,T < 1, but, in what follows we only need that Ds < 3 in this case.) 

If the tilt angle is made to be zero (£ = 0), we find that Ds is not given by 
(p3[), but by the following formula (which we derive in Sec. 7.2), 



D s < 3 for/i0 <h g + 1/t. 



(34) 



Since Ds is given by (|^) only if the tilt angle £ is precisely zero, we say that 
the untilted ellipsoid scattering system is atypical. As conjectured in ||, Ds 
from ( |33|) i s greater than or equal to Ds from (|34|). Note that Ds from the first 
line of (|33|) is larger that Ds from the first line of ( p4] ) by the factor h z /h^,. 
Although the transition of Ds from £ = to £ ^ is strictly discontinuous, 
there is also a continuous aspect: In numerically calculating the dimension of 
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a measure one typically plots ln/(e) versus ln(l/e), where /(e) = 5^^jln[l//Zj] 
and /ii is the measure of the cube in an e grid. One then estimates the 
dimension as the slope of a line fitted to small e values in such a plot. In the 
case of very small tilt, such a plot is expected to yield a slope given by (34) for 

e > and subsequently, for e < e*, to yield a slope given by (33). Here is a 
small tilt-dependent cross-over value, where e* — ► as £ — * 0. In such a case, 
the dimension, which is defined by the e — > limit, is given by (0). 



5 Numerical Computations for the Three-Dimensional 
Billiard Scatterer 

To verify that the untilted system is atypical we numerically calculated the 
box-counting dimension of SM for various values of the parameter rj_, then 
introduced a small tilt perturbation in the form of a — 27r/100 radian tilt of the 
ellipsoid about the x-axis and repeated the dimension calculations. The results 
are shown in Fig. ^. They confirm Eqs. ( |33| ) and (|3~4|). (We do not plot data 
points near the kink in the theoretical curve because of the numerical difficulty 
in obtaining reliable results in that parameter range. This is a result of slow 
converge near the kink. It is also found in the example of Sec. 3.3 and discussed 
in Appendix A.) 

We compare measured values of the box-counting dimension to the predicted 
values of the information dimension in Fig. ^. The box-counting dimension gives 
an upper bound on the information dimension, but often the values of the two 
dimensions are very close. In particular, we can compare the result for the tilted 
ellipsoid system to the 2D map of Sec. 3. In the regime h§ > 1/r, our system 
is similar to case (a) studied in Section 3.2 (note that, as discussed in Sec. 4, 
1/r is never smaller than h z ). Changing rj_ while leaving fixed changes r 
while h(f, changes only slightly. This is similar to varying Ai of the 2D map while 
keeping A2 fixed, i.e. varying r. Figure ||, shows (for a particular value of r) 
that the measured values of Do and D\ are numerically indistinguishable in the 
region corresponding to case (a), A < A a . 

For h^r < 1 {htf, < h z + 1/r) the information dimension of SM is predicted 
to be less than three for the tilted (untilted) ellipsoid system. This is analogous 
to cases (b) and (c) of Sec. 3.2; see Fig. |[ Since SM divides the four- 
dimensional phase space, its box-counting dimension cannot be less than three, 
and, analogous to the result of Sec. 3 and Fig. |], we expect that the box- 
counting dimension is three for h<f,T < 1 (h^ < h z + 1/t). This expectation is 
borne out by the numerical results, Fig. ^. 

The box-counting dimension of SM was computed using the uncertainty 
dimension method |], [Q). This method gives the box-counting dimension of 
the basin boundary which, as discussed above, coincides with SM. 

The uncertainty dimension method was carried out as follows: 

1. Choose a point, x, at random in a region of a 2D plane intersecting the 
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basin boundary (e.g., Fig. 6) and determine by iteration in which basin 
it lies. 



2. Determine in which basins the perturbed initial points x ± <5 lie (5 is some 
small vector). 

3. If the three points examined in (1) and (2) do not all lie in the same basin, 
then x is called "uncertain" . 

4. Repeat 1 to 3 for many points x randomly chosen in the 2D plane, and 
obtain the fraction of these that are uncertain. 

5. The fraction of points which is uncertain for a given 8, denoted /, scales 
like jL2| / ~ |<5| 2_do , where do is the box-counting dimension of the inter- 
section of SM with the 2D plane. The box-counting dimension of SM in 
the full 4D state space of the map is Do = 2 + do- (The dimension of a 
generic intersection of a 2D plane with a set having dimension Do i n & 4D 
space is c?o = 2 + Dq — 4, which gives Do = 2 + do.) Thus, plot In / versus 



In 



, fit a straight line to the plot, and estimate Dq as 4 minus the slope 



of this line. See Fig. @ for an example of such a plot. 



6 Structure of the Stable Manifold 
6.1 Untilted Ellipsoid 

Due to the symmetry to the untilted ellipsoid billiard, the chaotic saddle of 
the untilted ellipsoid system has a special geometry (i.e., it lies in z = v = 0), 
which, as we show, accounts for the dimension being lower than the predicted 
value for a typical system. Similarly, the symmetry induces a special geometry 
on the stable manifold (SM). Figure fj] shows a 3D slice of SM in the atypical 
system. The slice is at a fixed value of u>. The axes are z, v, and <j), but one 
could have chosen an arbitrary line through (4>,lu) as the third axes and seen 
a plot which was qualitatively the same. The stable manifold is organized into 
rays emanating from the 0-axis with oscillations along the (^-direction. The 
magnitude of the oscillations decreases to zero as the <^-axis is approached. 

To understand the structure of SM in more detail, assume that \z\ is small. 
Then, since \z\ — is invariant, we can approximate the dynamics by expanding 
to first order in z, 

z n+1 = DM z (0j n )z n , (35) 

$ n+1 = M^oJn). (36) 

Say (zsMi&Sm) is a point on SM. As this point is iterated we have that 
|z| — ► with increasing n. However, since ( |35| ) is linear in z n , for any constant 
a, and the initial condition, (olzsm,4>Sm), the subsequent orbit must also have 
\z\ — > 0. Consequently, if (zsm,4>sm) lies in SM, so does (azsM,4>SM)- Thus, 
for the system (^5|), (|36|), the stable manifold at any point <fi lies on a straight 
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line through the origin of the two-dimensional z-space. Put another way, in 
the approximation ([35]), ([36]), the stable manifold can be specified by an equa- 
tion giving the angle of z as a function of <fi. Thus, decomposing into polar 
coordinates (p, x), the stable manifold for \z\ — > approaches the form 



X = x(<t>)- 



(37) 



(This explains the structure seen in Fig. [?].) For \z\ finite the linearity of ( |35| ) is 
not exact, and we expect that the behavior of the stable manifold at constant (f> 
is not a straight line through the origin of the z plane. Rather, as \z\ becomes 
larger we expect (and numerically observe) the straight line for small \z\ to 
appear as a smooth curve through z = 0. 

Since for fixed <f> SM varies smoothly with increasing p = \ z\ , the dimension 
of SM is not affected by the approximation ( |35| ) and (36). That is, to find the 
dimension of SM, we can attempt to find it in the region of small \z\ where ([35]) 
and d36| ) are valid, and that determination will apply to the whole of SM. 

The task of analytically determining Dg is too hard for us to accomplish in 
a rigorous way for the system, (|35|), ([36]), applying to our billiard [a heuristic 
analysis yielding (|34|) is given in Sec. 7.2]. Thus, to make progress, we adopt 
a model system with the same structure as ([35]), (36). In particular, we wish 
to replace the two-dimensional billiard map (36) by a simpler map, M^ — > M^, 
that, like the original two-dimensional billiard map, is chaotic and describes a 
Hamiltonian system. For this purpose we choose the cat map, 



4> n+ l = M^n) = C<j> n modulo 1, 



(38) 



where C is the cat map matrix, 



C = 



2 1 
1 1 



Similarly, we replace DM z (0,(f>) in 
on </>, 

z n+1 = M z (z n ) = 



by a simple symplectic map depending 



A /(</>„) 
A- 1 



(39) 



where A > 1, and f(4>) is a smooth periodic function with period one in </> and 7r. 
[Note that vertical (i.e., parallel to z) line segments are uniformly expanded by 
the factor A, and thus, by the same argument as in Sec. 3.1, we have 1/r = In A 
and h z = 1/t.] 

Since only the angle of z is needed to specifiy the stable manifold, we intro- 
duce the variable v = z/v = tan%. We can then derive a map for v: From ([39]) 
we have v n+ iv n+ i = \v n v n + v n f(<fi n ) and v n+ \ = \~ 1 v„. Dividing the first 
equation by the second, and v n cancel and we obtain 



u n+ i = \ 2 v n + Xf '(</>). 



(40) 
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We now consider the dynamical system consisting of ( |38| ) and fl40|). Note that 
the system (|38|), ( f40| ) is a three-dimensional map, unlike the system (^), (|39|), 
which is a four -dimensional map. 

For (|3^), (|40|), the stable manifold is given by 

oo 

^ = ^(<?) = -A- 1 ^A- 2 7(C>o), (41) 

i=0 

where C is the cat map matrix. 

To verify that this is SM we note that points above SM, v > vs{4>) [below 
SM, v < v>s{4>)] are repelled toward v — > oo (2/ — > —00). Thus, on backwards 
iteration, points go toward SM. We take advantage of this behavior to deter- 
mine SM. Imagine that we iterate <f> forward n iterates to <f> n — C n <j) modulo 1, 
then choose a value of v n , and iterate it backwards using (^) . We find that the 
initial value of v at time zero giving the chosen value v n at time n is 



= K/A 2 ") - A -1 \~ 2l f{C l 



Keeping v n fixed and letting n — > +00, the value of vq approaches vs{<t>): given 
by ©. 

Results proven in ]l0| show that the box-counting dimension (the capacity) 
of the graph of the function v — vs{4>) given by (^l|) is 

D S = { 3-2^A, forA<5, 
s \ 2, for A > B, 1 ' 

where B — > 1 is the larger eigenvalue of the matirx C. The formula for 

A < B holds for almost all (with respect to Lebesgue measure) values of A. Since 
A > 1, the sum in converges absolutely, implying that vs{<j>) is a continuous 
function of (p. Thus, when the first result in Eq. ( |42"| ) applies (i.e., the surface 
is fractal with L>s > 2), the stable manifold is a continuous nondifferentiable 
surface. 



For example, evaluating (41) on the surface <j) = sii + where u + is the unit 
vector in the eigendirection of C corresponding to the expanding eigenvalue B, 
we have that v versus s is of the form 



i=Q 



Thus, the graph of v versus s for B > A has the form of Weierstra/3's famous 
example of a continuous, nowhere-differentiable curve. 



To obtain (42) in another way, we again consider the map (|38|), (|40|). We 
claim that 

can be regarded as a typical system, and that the formulae 
of Sec. 2 should apply to it. That is, in contrast to the existence of a symmetry 
for ( |39| ) [namely, z — > —z leaves ( |39| ) invariant], ( f40| ) has no special symmetry. 
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The Lyapunov exponent corresponding to (flO|) is h v = 2 In A. Note that for the 
system (p8|), ( p0| ) [and also for the system (|38|), (|39|)] there are no fluctuations 
in the finite time Lyapunov exponents and thus the decay time for the system 
(|38|), ( [40| ) is given by r^ 1 = h v . Noting that the Lyapunov exponents for the 
cat map are ± In B and applying (Q) to the three -dimensional map ([38|) , ([h^) , 
we immediately recover Eq. ([f2|). As discussed in Appendix B, this point of 
view can also be exploited for the original ellipsoid system [rather than just for 



the model system fl38| ) and (p9|)1 . 

Returning now to the full four-dimensional system, (|38|), (|39|), and noting 
that SM is smooth along the direction that we eliminated when we went from 
(j38|), ( |39|) to (|38|), (|4(i|), we have that the dimension D$ of the stable manifold of 
the invariant set (z = 0) for (|38j), ( |39] ) is £>s = + The Lyapunov exponents 
for z motion in the four-coordinate system are ±/i z = ± In A and ±/i0 = ± In B 
for (j) motion. In terms of the Lyapunov exponents, Dg is 

4-2£j, for h a /h+ < 1/2, 



£) J h*> 'i f - 1 > (43) 

\ 3, for h g /h+ > 1/2. V ; 

Since /i z = 1/r for (||) and (|9|) we see that, for h z jh^ < 1/2, Eq. (H|) is 
the same as Eq. (|3^). Also, the system (|3^), (|39| ) has no finite time Lyapunov 
exponent fluctations, and, thus, the information dimension and the box-counting 
dimensions are the same. Hence, D$ = 3 when SM is smooth ih, z jh^ > 1/2). 
This is analogous to the situation r = 1 and Eq. ( |2~i| ) of Sec. 3. [Th e lack 



of finite time Lyapunov exponent fluctuations for our model system fl38|), (|39|) 
is reflected in the fact that all periodic orbits in A have precisely the same 
Lyapunov exponents, namely ±lnA and ±ln£?.] 



6.2 Basin Boundary for a Map Modeling the Tilted Ellip- 
soid Billiard 

We now wish to investigate the structure of the stable manifold when we give 
the ellipsoid a small tilt. Again, we adopt the approach of Sec. 4.2: we obtain a 
rigorous result by utilizing a simpler map model that preserves the basic features 
of the tilted ellipsoid case. Here we again use ([38]) but we now modify (|3^) to 
incorporate the main effect of a small tilt. The effect of this modification is to 
destroy the invariance of z = 0. Thinking of the first non-zero term in a power 
series expansion for small |z|, this invariance results because the first expansion 
term is linear in z\ i.e., the z- independent term in the expansion is exactly zero. 
When there is tilt this is not so. Thus we replace ([39]) by 













+ 


' fz{$n) ' 


V n +l 




A- 1 








fv(<t>n) 



(44) 



The simplest version of (Q) which still has the essential breaking of z — > — z 
symmetry is the case where / = f v = 0. Because setting / = f v = greatly 
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Zn+1 




'A 








fz{<f>n) 


Vn+1 




A- 1 _ 






+ 






simplifies the analysis, we consider this case in what follows (we do not expect 
our conclusion to change if /, /„ ^ 0). Thus, we have 



(45) 



The problem of finding the stable manifold for the invariant set of the map, ( p8[ ) 
and ([l5|), is now the same as for the previously considered case of Eqs. ( |38| ) and 
( ^0[ ) [compare the equation z n+ i = Xz n + f z {4>n) with Thus, making use 

of this equivalence we can immediately write down the equation for the stable 
manifold in the four-dimensional state space (z,v,<fi,(jj) as 



DO 

■E 

i=0 



z = -vY\- i f z (C i d>), 



(46) 



which is obtained from ( |4l| ) using the replacements v — ► z/u, A/ — > / z , and 
A 2 — > A. The rigorous results of Ref. jl0| again show that this is a continuous, 
nowhere-differentiable surface for almost all A in 1 < A < B, and, furthermore, 
when this is so (In A/ In .B = h z /h<f, < 1) we have 



D s 



hd 



(47) 



Also, Ds = 3 when hz/h^, > 1. Since h z = 1/r, this is the same as Eq. (33). 



7 Derivation of Dimension Formulae 

7.1 Ds for Typical Systems 

While the dimension formulae presented in Sec. 2 and derived in were found 
to be valid for the two-dimensional map of Sec. 3 and the ellipsoid scattercr of 
Sees. 4-6, the derivation given in Q does not apply to these systems. In the 
derivation of M, a higher-dimensional generalization of horseshoe- like dynam- 
ics is assumed. When this type of dynamics is responsible for a fractal basin 
boundary, the boundary looks, locally, like a Cantor set of smooth surfaces. The 
examples of this paper do not possess horseshoe-like dynamics. Their dynamics 
produces basin boundaries which consist of a single continuous surface; see Figs. 
0, ^|, and 0. We proceed to derive dimension formulae using a more general ar- 
gument not restricted to an assumed particular type of boundary. The resulting 
formulae are identical to those derived in ^ and stated in Sec. 2. 

Let R be the portion of the state space which contains all points within e 
of a nonattracting, ergodic, invariant set, A of an M-dimensional map, P. Let 
sm be the portion of the stable manifold of A which is contained in R. The 
points in R are within e of sm since A is a subset of sm. If we sprinkle a large 
number, A^o, of orbit initial conditions in R, then the number of orbits left in R 
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after n » 1 iterates is assumed to scale like (||) N ti /Nq ~ e ™/ r . Let the map, 
P, have Lyapunov exponents 

hy > h^ J _ 1 >--->hf>0> -hi >■■■> -hs_x > -hg, 

where U + S = M. 

We will use the box-counting definition of dimension 

N(e) ~ e~ D , (48) 

where N(e) is the minimum number of M-dimensional hypercubes ("boxes") 
needed to cover sm and D is its box-counting dimension. 

We wish to develop a covering of sm using small boxes and determine how 
the number of boxes in that covering scales as the size of the boxes is decreased. 
We will look at how the linearized system dynamics distorts a typical small 
box. This will help us determine how the Lyapunov exponents are related to 
the dimension of sm. Since we assume a smooth map, the dimension of the 
stable manifold of the map is equal to that of sm. 

Cover sm with boxes which are of length e on each of their M sides. Call 
this set of boxes C. The number of boxes in C is denoted Nq. Iterate each 
box forward n steps, where n is large, but not so large that the linearized 
dynamics do not apply to the boxes. Call the set of iterated boxes P n (C). A 
typical iterated box in P n (C) is a distorted Af-dimensional parallelopiped and 
has dimensions 

ee nh i x • • • x ee nh l x ee~ nh ^ x • • • x ee~ nh s . 

Each parallelopiped intersects sm since its preimage (a box) did. (The set 
P n (C) does not, however, cover all of sm since points on sm move closer to A 
on forward iteration.) 

To construct a refined covering of sm we begin by covering each parallelop- 
iped of P n {C) with slabs of size 

u factors 

e x ■ • ■ x e xee~"' l i x ■ • • x ee~ nh s . 

There are roughly exp[n(hy + • • • + h^)] such slabs. Only ~ e~ n l T of these slab 

are within e of sm. Let C denote the set of N' ~ No exp[n(h^} H h h{ — 1/r)] 

slabs needed to cover the part of sm lying in P n (C). Let us define H — + 
■ ■ ■ + ^ - 1/r so that N' - Aoe' 1 ^. 

Iterate each of the N' slabs backward n steps. We now have the set P~ n (C). 
It contains N' parallelopipeds of size 

+ + S factors 

ee ~ nh u x • • • x ee~ nhl X e X • • • X e . 

The set P~~ n (C) forms a covering of sm. To calculate the dimension (which 
is defined in terms of boxes) we cover the N' parallelopipeds in P~ n (C r ) with 
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boxes which are e 3 = te~ nh i+ x on each side. (We will choose the value of the 
index j below.) The number of boxes needed to cover sm, for boxes of size tj, 
scales as 

„.// , T / ee 3 ee ee i / 1 
N ~ N' x x • • • x > 



c 3 C J c 3 \ c 3 . 

To cover the slabs with these boxes we need a factor of 1 boxes along each 
direction of a slab which is shorter than the edge length of a box and a factor of 
ee~ nh k /ej boxes along each direction (here, the kfi 1 direction) which is longer 
than the edge length of a box. In terms of No, 

N" ~ N exp{n[(S + j)h j+1 + H - hf - hf hf - 1/r]}. 

To compute the dimension of sm, we compare N(e) = N to N(ej) = N" 
using (ktq). This gives N(ej)/N(e) ~ {e/ej) D ~ exp(n£>/i+ +1 ) which yields the 
following j-dependent dimension esitmate, 

H -hf + h+ + ■■ -hf 
D(j)=S + j+ (49) 

Our definition of box-counting dimension, (|4^) , requires us to find the minimum 
number of boxes needed to cover the set. Since we are certain that the set is 
covered, ( |4^ ) yields an upper bound for the dimension for any j. Thus, to find 
the best estimate of those given by ([49|), we minimize D(j) over the index j. 
Comparing D(j) to D(j + 1) yields the condition 

hf H h hfj + hj > H > hf H \-hj, 

where J is the best choice for j (i.e., the choice giving the minimum upper 
bound). The conjecture is that this minimum upper bound D(J) from (E^) is 
the actual dimension of sm for typical systems. D( J) is the same as the result 
Eq. (|^) presented in Sec. 2 for D$- The derivation of the dimension formula 
for the unstable manifold, Djj, is similar to that just presented for the stable 
manifold. 

[The derivation just presented gives the information dimension of the sta- 
ble manifold, not the box-counting dimension. We were considering the sizes 
of typical boxes in the system and covered only those. This leaves boxes that 
have atypical stretching rates for large n (a box containing a periodic point, for 
example, will, in general, not stretch at the rates given by the Lyapunov expo- 
nents) unaccounted for. What we have actually computed is the box-counting 
dimension of most of the measure, which is the information dimension ]Tl[ | , of 
the stable manifold. See for discussion of this point.] 

As derived, D(J) gives an upper bound on the dimension. We saw in the 
ellipsoid example of Sees. 4-6 that the z — > —z symmetry of the system lead to 
a special geometry for its stable manifold (Fig. 0) and the formula just derived 
did not apply (although it was an uppder bound). It is the conjecture of Q 
that this formula gives not the upper bound, but the exact dimension of the 
stable manifold for typical systems. This is supported by our results for the 
tilted ellipsoid example. 
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7.2 D s for the Untilted (Atypical) Case 

The derivation above for typical systems gives the wrong dimension formula 
for the stable manifold of the (atypical) untilted ellipsoid system. This typical- 
system derivation consistently overcounts the boxes needed to cover the atypical 
stable manifold of the untilted ellipsoid system. Sec. 6.1 showed that the 
stable manifold had a ray-like structure (cf. Fig. |?|) induced by the z-direction 
reflection symmetry of the system. If we account for this structure we can 
modify the typical-system derivation so that we do not overcount boxes. The 
modified derivation produces the correct formula. 

Recall that our untilted ellipsoid map (which we call P here) has Lyapunov 
exponents ±/i z and ±/i0 and decay time r. We consider the case 

h<j, > h z + t > 

Let us cover the region of state space which is within e/2 of A with No 
boxes with edge lengths e x e x e x e. Call this set C. The map, P, is linear 
in z (in particular, the z portion is of the form z n +i = DM((p)z n ) so that the 
graph of SM in z space for a given value of (f> is a straight line through z — 0, 
i.e. z = vg(4>). We denote by sm the portion of the stable manifold which is 
contained in C. Since SM contains A, each box in C intersects sm. 

Iterate these boxes forward n ^> 1 steps. They become a set distorted of 
parallelopipeds which we call P n (C). Each parallelopiped has dimensions 

ee nh„ x ee nh z x ee -nhz x gg-nft* 

and intersects sm. We cover P n (C) by N exp[n(h z + h^)] parallelopipeds with 
dimensions 

e x e x ee- nh * x ee""' 1 *, 

and discard all of these parallelopipeds which do not intersect sm. There are 
N' ~ Ao exp[n(h z + h<p — 1/t)] parallelopipeds remaining which cover sm. The 
portions of these paralleopipeds which are in e > \z\ > ee~ nh * do not contain 
sm. Suppose some portion of sm did fall in e > \z\ > ee~ nhz . Since P is 
linear in z, this outlying portion of sm would, upon n reverse iterations, map to 
the region ee nhz > \z\ > e, which contradicts the definition of sm given above 
(i.e. sm is within e/2 of A). Therefore we can discard the portion of each 
parallelopiped which lies is e > \z\ > ee~ nh * . These N' parallelopipeds now 
have dimensions 

e x ee- nh * x ee- nhz x ee~"' 1 *, 

and are denoted C . We now iterate the parallelopipeds in C backward n times 
and call the resulting set P~ n (C). This set contains N 1 parallelopipeds with 
dimensions 

ee"™' 1 * x ee- 2 "' 1 - x e x e. 

We cover P~ n (C) (which covers sm) with hypercubes which have as their edge 
length Ee~ nh * . The number of hypercubes needed is A^" ~ Noexp[n(h z + — 
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1/t — 2h z + 3h$)]. The (information) dimension is again found by comparing 
the number of e-sized hypercubes needed to cover sm to the number of ee~ nh * 
sized hypercubes needed to cover sm, 



N" fee- nh " 



— D 



which yields 

D = 4 



h z + 1/t 



h<j, 

This expression is between three and four for (h z + l/r)/h^ < 1. 

8 Discussions and Conclusions 

In summary, this paper has tested and illustrated formulae giving the fractal 
dimension of nonattracting chaotic sets in terms of their Lyapunov exponents 
and decay time. We have done this using two examples, one, a noninvcrtiblc, 
two-dimensional map with two expanding (positive) Lyapunov exponents, and 
the other, a chaotic scattering billiard in three spatial dimensions. The first 
example is particularly useful for understanding the transient measure. The 
second example provides a striking illustration of the issue of atypicality, and 
also has the benefit of potential physical realization (see below) . Another point 
is that our second example provides the first known case of a chaotic scatterer 
with a fractal basin boundary that is a continuous, nowherc-diffcrentiable sur- 
face (such structures have been previously discussed for basins of attraction of 
dissipative systems but not in chaotic scattering). Finally, in Sec. 7 we have 
provided arguments that are more general than those previously given Q for 
the dimension formulae of Sec. 2. 

To conclude, we now offer some discussion on the possible experimental real- 
ization of our billiard example. Our point here is that such billiard systems offer 
a particularly nice way of experimentally realizing and studying basin bound- 
aries. In the past there have been very few experiments that have attempted to 
study basin boundaries. The problem is the experimental difficulty of carrying 
out the technique used in numerical work: namely, choosing a grid of initial 
conditions and seeing where each orbit from each initial condition goes. To do 
this experimentally, one would have to precisely prepare an initial condition, 
run the experiment to see where the orbit with that initial condition goes, and 
do this for each initial condition. The difficulties with this approach are that it 
is often not possible to prepare initial conditions either on a fine enough grid or 
sufficiently precisely to observe small scale basin boundary structures and that 
running the experiment many times can be time consuming. (In addition, exper- 
imental parameters may drift over the course of the multiple experimental runs.) 



Nevertheless, in one case |13 this program was successfully carried out. In that 
paper an electrical circuit was used as the experimental chaotic system, and the 
authors were able to experimentally demonstrate a riddled basin of attraction 
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(a type of basin for which final outcomes are particularly difficult to predict; 
see Refs. Q, In another work, Cusumano and Kimble [^6j have devised 

a new experimental procedure for studying basin boundaries. This procedure 
has advantages over the straightforward method of preparing initial conditions 
on a grid. Using this method, Cusumano and Kimble have successfully mapped 
out basin boundaries for experimental dissipative mechanical systems. To our 
knowledge, these two works (and a billiard we discuss below) are so far the only 
experimental investigations of fractal basin boundaries. This is in contrast to 
the large number of papers that have experimentally realized fractal structure 
in chaotic attractors. This disparity in the situations of attracting and nonat- 
tracting chaotic sets seems to be largely due to the disparity in the ease of 
experimental realization for the two cases. 

We believe that billiard systems offer a convenient avenue for experimental 
investigation of the various general types of basin boundaries. In particular, we 
suggest an optical realization of chaotic scattering billiards in which the billiard 
surface is made to be optically reflecting, and light rays play the role of orbits. 
For the scatterer of Sec. 4-6, for example, imagine that we have a tube with 
cross section as shown in Fig. ^(b) and mirrored inner walls. The length of the 
tube is L > 2r±, where r± is the major radius of an ellipsoid with a mirrored 
surface. Imagine that the tube is oriented vertically with its bottom, open end 
placed on a red surface, and that the ellipsoid is suspended in the tube. In this 
configuration, an observer looking in the top of the tube sees multiple reflections 
of the red surface that is at the bottom of the tube. Since rays whose directions 
are reversed retrace the same path, we can think of orbits from initial conditions 
starting on the retina of the observer's eye passing through his pupil, bouncing 
around in the scatterer, and then exiting either through the bottom (red) or the 
top (not red). Thus, the observed boundary of the red region is precisely the 
basin boundary that we have been discussing in Sees. 4-6. Note that to observe 
that boundary we did not have to prepare many initial conditions and repeat 
the experiment many times. Rather a single image giving a global picture of the 
basin can be formed (replacing the observer's eye by a camera). We have already 
applied this approach to a chaotic scatterer formed from four mirrored spheres 
to experimentally study a different general type of basin boundary structure 
known as a Wada boundary (l?]] (a Wada boundary is a boundary separating 
three or more regions such that every boundary point is a boundary point for 
all regions). In the case of our ellipsoid billiard, such an experiment would 
represent the first experimental realization of a fractal basin boundary with 
the structure of a continuous, nowhere-differentiable surface, one of the small 
number of typical possible basin boundary types . 

We thank J. A. Yorke and B.R. Hunt for discussion. This work was supported 
by ONR and DOE. 
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Appendix A: Slow Convergence of the Dimension 
Near the Transition Between Fractal and Non- 
fractal Behavior 

It is difficult to numerically measure the box-counting dimension of the chaotic 
repellors studied in Sec. 3 and Sees. 4-6 for some range of the system parameter. 
When the parameter is set near the value at which the repellor makes the 
transition from a fractal to a non-fractal function the number of boxes needed 
to cover the set scales as 

»«>~^ (AD 



with e (the width of a box) down to very small e. In ( Al), Dq = 1 in the case of 



the two-dimensional map and Dq = 3 in the case of the ellipsoid system. Due 



to the logarithmic term in (Al) accurate numerical determination of Dq can 



be very demanding. In the following we demonstrate that (Al) holds near the 
transition point between a fractal and a nonfractal. In particular, we consider 
[-Do = 1 in Eq. (Al)] in the case where Ai = A 2 = A. Using the same 



reasoning as for the derivation of Eq. ( pfl| ) we obtain for the invariant set the 



curve 



*) = -^EV lsil >( 2A )- (A2) 

i=0 

This result is also derived in Ref . [j]] . The box-counting dimension of the graph 
of y{x) versus x is @ D Q = 2 - [In A/ In 2] for A < 2 and Dq = 1 for A > 2. 
Thus, the graph of y(x) is a fractal for A < 2 and is non-fractal for A > 2. We 
are interested in the transition region, A = 2, for which we wish to show that 



(Al) applies. 



To calculate the box-counting dimension of (A2) we need to know how the 
number of e-width boxes needed to cover the set scales with e. If we let the 
sum in ( |A2j ) run from i = to i = n, the resulting approximation to y(x) is a 
smooth, finite length curve: 

n 

Vn{x) =~Y\ A~ l sin(27r2^). (A3) 

i=Q 

As n increases, the ratio of the height to the width of the oscillations added 
by each term in the sum increases. For the purpose of estimating N(e) for 
some fixed value of e it suffices to consider ( [A3] ) with n sufficiently large. A 
practical size for the boxes in a covering is one which scales as the width of the 



smallest oscillations of (A2): e n ~ 2 71 . This wa y we can be sure to resolve the 



smallest details of the approximation curve (A3) and be sure that higher order 
approximations will not significantly alter our count of the number of boxes in 
the covering. (Oscillations introduced by the n + l s ^ term, for example, will be 
roughly half the width of a box of size e„ ~ 2~™, or ~ 2~(" +1 ) .) 
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The length of the curve ( |A3| ) is 

£ n = J dx^iTWM 2 = (V^+EM?) > 

where (•) denotes an average over x. For large n and A < 2, is typically 

large and i n (\y'(x)\). 

To determine £ n we examine 

" / 2 \ 1 

i=0 ^ ' 

For j sufficiently large, the largest term in this sum dominates and y' n (x) « 
(|) n . In this case, £ n ~ (f)" and N(e n ) — |a. Taking e = 2 _n , or n = 
— Ine/ln2 we have that N(e) ~ e -(2-iaA/in2) < rpj^ sca j m g j g a p 0wer j aw j n e . 
(It gives a box-counting dimension of Dq = 2 — In A/ In 2.) 

For A < 2, but A w 2, the last term will not dominate the sum, but the 
last several, up to n x , terms will dominate. To estimate n x , set (jj x equal to 
a constant factor if > 1 where K — 1 is order one (e.g., X = 2). This gives 
n x ~ (2 — A) -1 ^> 1. If n < n x , then all the coefficients (2/A) 1 in the sum 
y' n (x) are of the same order, y' n (x) ~ X)"=o cos(27r2*x). For a given value of x 
and different values of i ^> 1 the arguments of the cosine terms in the sum are 
very different and vary rapidly with x. Thus, we regard the terms in the sum 
as random. In this case y' n ~ \/n, and, thus, i n ~ \fn. Again, we say that 
N(e n ) ~ j^, but now N(e n ) — \fnje n — V- hie„/e„, which is (|A]]). If n > n x , 

thenK(x)|~ni /2 (f) n . 

For the two-dimensional example studied numerically in Sec. 3.3, the numer- 
ical results are shown in Fig. |^. Due to the effect discussed in this appendix we 
have not plotted data near the kink at Ai = A&. The results for the ellipsoid sys- 
tem are discussed in Sec. 5 and shown in Fig. |6|. A typical two-dimensional slice 
of the basin boundary for this system is a continuous, nowhere-diffferentiable 
curve similar to the chaotic saddle of the two-dimensional system. We thus 
have the same difficulty in measuring the dimension near h^r = 1, Fig. ||(a) 
and = h z +T- 1 Fig. |(b). 



Appendix B: Dynamics of the Orientation of z 

Rewriting ([35]) as, 

Z n+1 = mii(<f>n)Zn + mi2(<f>n)v (Bl) 
Vn+1 = TO2l(^»)z n + m 2 2(^n)«n, (B2) 

eliminating z in favor of v — z/v, and dividing the first of the above equations 
by the second, v„+i and v n cancel to yield an evolution equation for v, 

mn(<f>n)Vn + Wl2(0n) /T30 s 
v n+l = ~J- -j—. (ad) 

m 2 l((f> n )Vn + m 22 ((Pn) 
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Equations (J3q) and (B3) now constitute a 3D map. Our transformation from 



(Bl), (B2) eliminates the symmetry of the original system, and it is reasonable 
to now assume that the information dimension of the natural transient measure 
for the 3D map (p6|) and (BS) is given by the result for a typical system. Thus, 



we have that the information dimension of the stable manifold of the invariant 
set for the 3D map is 



£»s = 3-(^f)- 1 ! for/i f>l. 



(B4) 



Where the superscribed circumflex denotes quantities calculated for the nat- 
ural transient measure of the 3D map, ( |36| ) and (BS). Since the information 
dimension D$ is a lower bound on the box-counting dimension for the map 
( |36|) and ( |B3[ ) , and since the box-counting dimension of SM for the original 4D 
map is one plus that for the 3D map (this follows from the structure shown in 
Fig. [?]), we have that £>s + 1 [with L>s given by (B4)] is a lower bound for the 
box-counting dimension of SM for the 4D map. We emphasize, however, that 
£>s + 1 is different from D$ for the 4D system. This is because of the differ- 
ence in the natural transient measures for the 3D and 4D maps. (In contrast, 
the box-counting dimensions do not depend on the measures.) In particular, 
as shown below, as compared to the 4D map natural transient measure, the 
natural transient measure for the 3D map more strongly weights <j> values which 
experience slower repulsion from the invariant set A. 

To see this we consider the Lyapunov exponent corresponding to the evolu- 



tion of a differential perturbation in v which we denote 8v. Differentiating (B3) 
we have 

5v n+x /5v n = d($ n ){m 2 iv n + m 22 )~ 2 , (B5) 



where d{cf> n ) the determinant of (Bl), (|B2|) , d(cf> n ) = miiTO 2 2 — ?«i2 TO 2i = 1- 
Noting that v n = z n /v n and comparing ( |B5[) with (B2), we have 



8v n+1 /8v = d(4> n ){v n /v n+1 ) 2 . 

For points on SM , v n decreases exponentially with time (i.e. 
invariant set A). Hence, 



5v n /5is = Q(v n /v ) 



exp[2nh z (<f>o,ri)], 



(B6) 

it approaches the 



(B7) 



where Q = d((f>o)d(<fii)...d(<fi n -i) and h z (<p,n) with = (f>Q is the expanding z, 
finite-time Lyapunov exponent for the 4D system. The estimate, exp(2n/i z ), in 
(B7) follows since, due to the Hamiltonian nature of the problem, Q 1 ^ 71 1 as 
n — > 00. Thus, for a given value of <p 1 in the 3D system, the expansion away 
from the invariant set is at the rate 2h z ((p,n) rather than at the rate h z ((j),n), 
applying to the 4D system. For the example of ( |39| ) there are no finite time 
Lyapunov exponent fluctuations [h z ((f>, n) = In A independent of n and <fi], so 
that the 3D and 4D natural transient measures are the same, both are uniform 
in <fr. 

In this case, 

f = t/2 = (In A)" 1 . (B8) 
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More generally, the finite time exponents, h z (cf>,n), fluctuate as <f> varies. The 
faster expansion rate away from A (namely 2h z ) in the 3D case means that <fr 
values with slower expansion rates are more strongly weighted. Thus, for a case 
such as our ellipsoid problem we expect that 

f > t/2% v < 2h z , 

where h v is the infinite time Lyapuonv exponent for the 3D map measure and 
for perturbations in v. 
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Figure 1: Basins of 2D map model, Eqs. @-@ for Ai = 1.1, r = 3. A 
500x500 grid of initial conditions was chosen in < x < 1, —0.3 < y < .2. Each 
was iterated forward by the 2D map until the orbit fell in one of the regions 
\y\ > 5. If the orbit fell in y > 5 (y < —5) a white (black) point was plotted on 
the initial condition grid shown above. The dimension of the boundary between 
black and white (the basin boundary) is D « 1.28. 
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Figure 2: Comparison of prediction and numerical dimension calculations for 
the 2D model map, with r = A2/A1 = 3 fixed. The filled dots show 

calculations of the box-counting dimension, and the boxes show information 
dimension. The solid line shows the predicted information dimension as given 
by Eqs. 
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h* 1/t h h 

Figure 3: G(h) versus h. The long-dashed line has slope —1 and is tangent to 
the graph of G(h) at the point h = h*. 
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Figure 4: The scattering system in the (atypical) case of vertical orientation of 
the ellipsoid, (a) A hard ellipsoid is placed inside a hard tube with cross-section 
as shown in (b). (The circle in (b) is the cross section of the ellipsoid.) R = 25, 
d = 10 and the radius of the ellipsoid at z — is — 5. 
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Figure 5: Basins for z — > oo (white) and z — > — oo (black). r± = 7, y = 5.1, 
4^ = 0, Us = .1, and Wy is given by the condition \v\ = 1 for (a) the untilted 
case and (b) a tilt of . 
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Figure 6: Comparison of dimension formulae, Eq. (33) and (|34|), with numerical 
estimates. The lines are the graphs of the equations and the dots show the 
numerical estimates of the box-counting dimension, (a) The tilted ellipsoid 
(typical) case, (b) The untilted ellipsoid (atypical) case. 
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Figure 7: Stable mainfold (SM) of the chaotic saddle (A); r± = 15, y = 5.1, 
v x = 0, and v y is given by the condition \v\ = 1. The t> 2 = 0, z = line is 
A. Notice the organization of SM into rays emanating from A. This special 
structure is atypical and is used to derive the atypical dimension formula for 
SM, Eq. (||). 
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Figure 8: Uncertainty dimension measurement for the untilted ellipsoid system 
with rj_ = 100. Plotted is the fraction of a random sample of points which are 



-uncertain" for various values of 



. The line running through the dots is a 



linear fit with a slope of w .07, yielding a dimension of do ~ 2 — .07 = 1.93. 
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